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We consider a two-dimensional electron gas in the presence of Rashba spin-orbit coupling, and 
study the effects of magnetic s-wave and long-range non-magnetic impurities on the spin-charge 
dynamics of the system. We focus on voltage induced spin polarizations and their relation to spin 
Hall currents. Our results are obtained using the quasiclassical Green function technique, and hold 
in the full range of the disorder parameter appr. 
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In the field of spintronics, much attention has recently 
been paid to spin-orbit related phenomena in semicon- 
ductors. One such phenomenon is the spin Hall effect, 
i.e. a spin current flowing perpendicular to an applied 
electric field^'^'^'^. It is now well known that for linear-in- 
momentum spin-orbit couplings like the Rashba or Dres- 
selhaus ones the spin Hall current vanishes exactly in 
the bulk of a disordered two-dimensional electron gas 
(2DEG)^'^'^'^. This can be understood by looking at 
the peculiar form of the continuity equations for the 
spin, as derived from its equations of motion in oper- 
ator form^iiSiiii. For a magnetically disordered 2DEG 
things are however different, and a non- vanishing spin 
Hall conductivity is found —li^ii^. Once more, a look at 
the continuity equations provides a clear and simple ex- 
planation of the effect^: a new term, whose appearance 
is due to magnetic impurities, directly relates in-plane 
spin polarizations, induced by the electric field, to spin 
currents. As the former, which have been the object of 
both theoretical and experimental studie o^^'^^i^^'^^i^^'^° , 
are influenced by the type of non-magnetic scatterers con- 
sidered, we forgo the simplified assumption that these be 
s-wave, and take into account the full angle dependence 
of the scattering potential. Besides going beyond what is 
currently found in the literature, where, in the presence 
of magnetic impurities, the non-magnetic disorder is ei- 
ther neglected or purely s-wave, our approach also shows 
the interplay between polarizations and spin currents in 
a 2DEG^. We note that in the correct limits our results 



agree with what is found in Ref. . On the other hand 
a discrepancy with Ref. [13 arises. 

For the calculations we rely on the Eilenberger equa- 
tion for the quasiclassical Green function in the presence 
of spin-orbit coupling^^ . The spin-orbit energy is taken to 
be small compared to the Fermi energy, i.e. app ^ £f 
- or equivalently a <^ vp - and the standard metallic 
regime condition 1/r ^ ep is also assumed. Here a is 
the spin-orbit coupling constant, pp{vp) the Fermi mo- 
mentum (velocity) in the absence of such coupling, and 
r the elastic quasiparticle lifetime due to non-magnetic 
scatterers. Our results hold for a wide range of values 



of the dimensionless parameter appT, since this is not 
restricted by the above assumptions. Contributions of 
order {a/vp)"^ are neglected throughout. We focus on 
intrinsic effects in the Rashba model; extrinsic ones^^, 
Dresselhaus terms^i and hole gases^^ are not taken into 
account. Finally, weak localization corrections, which 
could in principle play an important rolei^, are beyond 
the scope of our present work. 

The Hamiltonian of the 2DEG, confined to the x-y 
plane, reads 



H=^-h-cT + V{^), 
2m 



(1) 



with b — aez x p the Rashba internal field, cr the vector 
of Pauli matrices, and F(x) ~ Vnm{^) + Kn(x) the disor- 
der potential due to randomly distributed impuritiesS^. 
Non-magnetic scatterers give rise to Vnm(x) 



Km(x) =^ ;7(x-R,), 



while Vin(x) describes magnetic s-wave disorder 



Kn(x) 



B • cr(5(x - R,). 



(2) 



(3) 



Both potentials are treated in the Born approximation, 
and the standard averaging technique is applied. 

To begin with, we look at the continuity equation for 
the Sy spin polarizationi^i^i 



dtSy -I- 3x • 



3t, 



(4) 



where the second term on the r.h.s. is due to magnetic 
impurities. Here Tgf is the spin-flip time which stems 
from the potential ([3|) [cf. Eq. (fT3|) ]. Under stationary 
and uniform conditions the above equation implies a van- 
ishing spin current - hence a vanishing spin Hall conduc- 
tivity - unless magnetic disorder is also present, in which 
case instead 



3maT. 
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ization is given by Sy 



'Ntia\e\£TE, with 
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FIG. 1: (Color online) (a), (b) - The Fermi surface shift, 
Sp = \e\£T, due to an applied electric field along the x- 
direction. The white arrows show the direction of the internal 
field b. (c), (d) - Shifted bands and spin polarization in sta- 
tionary conditions, (c) - Asymmetric shift of the two bands 
when angle dependent scattering is present. The long dark 
(blue) arrows show the contributions to the spin polarization 
arising from a sector d(p of phase space, (d) - When magnetic 
disorder is turned on, additional contributions orthogonal to 
the internal field b appear, here shown by the short inward 
and outward pointing (blue) arrows. Out-of-plane contribu- 
tions are also present, but for the sake of simplicity not shown. 



Since the out-of-plane polarized spin current is related to 
the in-plane spin polarization, we now use simple phys- 
ical arguments to explain how the latter is generated 
by an applied voltagei^'^®. Since the Fermi surface is 
shifted by an amount proportional to the applied elec- 
tric field (say along the x-direction) , as shown in Fig. [T] 
(a),(b), there will be more occupied states with spin up 
- along y - than with spin down. In the case of short- 
range disorder, the total in-plane polarization can be es- 
timated to be proportional to the density of states times 
the shift in momentum, Sy ~ NSp N\e\£T. Since 
in the present situation we are dealing with the two 
Fermi surfaces corresponding to the two helicity bands 
e± = /2mzLap, obtained from the Hamiltonian ([T]), one 
expects Sy ~ {N+—N^)Sp, where, for the Rashba interac- 
tion, one has N± = No{l =F a/vp), Nq = m/27r. Explicit 
calculations agree with this simple picture and lead to the 
result due to Edelstein^'"', Sy = —NQa\e\£T. When long- 
range disorder is considered, a reasonable guess could be 
to substitute for r the transport time Ttr 



Ttr, —= f deW{9)il-cos{9)), 

Ttr J 



(6) 



W{9) being the angle-dependent scattering probability, 
so that Sy = —NQa\e\£Ttr. This was proposed in [29| . 
however the picture is too simplistic, and therefore the 
guess is wrong. As discussed in [2l|, the proper Sy polar- 



TE, — ^ [ dew {9) {l-cos{20)). 
te J 



(7) 



This particular time te, where stands for Edelstein, 
arises from the asymmetric shift of the two Fermi sur- 
faces, as depicted in Fig.[T](c), due to different transport 
times in the two bands. It shows that contributions from 
both forward [9 = 0) and backward {9 = tt) scatter- 
ing are suppressed. The next step is to consider what 
happens when magnetic impurities are included. Rely- 
ing once again on the simple picture of the shifted Fermi 
surface, one could argue that these have a rather small 
impact on the spin polarization, since the spin-flip scat- 
tering time usually makes a small contribution to the 
total transport time. However, even when this is the 
case, magnetic disorder does not simply modify the total 
transport time, but has an additional non-trivial effect. 
In its presence the spins do not align themselves along the 
internal b field, since they acquire non- vanishing compo- 
nents in the plane orthogonal to it - see Fig. [T] (d). It 
is these components who give rise to a finite spin Hall 
conductivity. In this respect, magnetic disorder has an 
effect similar to that of an in-plane magnetic field: it af- 
fects the spin quantization axis and tilts the spins out of 
their expected stationary direction. We now make these 
arguments quantitative. 

The starting point is the Eilenberger equation^^ , which 
we write explicitly for a homogeneous Rashba 2DEG in 
linear response to a constant and homogeneous applied 
electric field 



1 K 



(8) 



The quasiclassical Green function {g = gt^t^ii^'i^)) is de- 
fined as (^ — p^ /2m — fi) 

5 = ^y"d^G'tit2(p,x), G = ^ G-^ ) ' 

where Gt^t2(p,x) is the Wigner representation of the 
Green function, which has a matrix structure in both 
Keldysh (denoted by the check symbol) and spin space. 
Eq. ([5]) is the equation of motion for the Keldysh com- 
ponent - the one related to physical observables - iden- 
tified by the superscript "if" , which will be from now 
on implicitly assumed and thus dropped. Moreover, 
g^q = tanh(e/2T)(5f^ - 5^), where g^^ = -g^^ = 
1 — i9jb • <T, indicates the equilibrium - no electric field 
- function^^. All objects are evaluated at the Fermi sur- 
face in the absence of spin-orbit coupling, while ip is the 
angle defined by the momentum, p — p(cos ip, sin ip), and 
= { — sin ip, cos ip). From Eqs. ([2]) and ([3]) one obtains 
the self-energy contributions 



Snm(p) = nnn,;^ |C/(p - p')l'G(p') 



(10) 



3 



and 



3 



(11) 



1=1 p 



where and denote the concentrations of non- 
magnetic and magnetic impurities, respectively. In order 
to consider long-range non-magnetic disorder, we first 
expand the non-magnetic scattering kernel in spherical 
harmonics of the scattering angle and neglect its depen- 
dence on the modulus of p and p' 



1 



(l + 2i^iCos(^-^') 
2'kNqt 

2i^2COs(2(^ - 2ip') + ...) 



(12) 



with T the non-magnetic contribution to the elastic life- 
time. Then we write the magnetic scattering kernel in 
terms of the spin-flip time r,/ 



1 



2'KNnT. 



(13) 



The complete disorder self-energy can then be written 
separating its s-wave and higher harmonics contributions 



' (a) - ^{Kg) (14) 



&Tsf ^ 2r 2r 

where (...) = J Atf/2Ti.... 

The connection between g and the physical observables 
is made by integrating over the energy e, which is the 
Fourier conjugate variable of the time difference ti — ^2- 
For instance, the spin density is given by the angular 
average of the Keldysh componenlj^ 



s = s 



^ j de(Tr(<T5)). 



(15) 



In order to solve Eq. (I8|), it is convenient to turn it into 
matrix form, writing g as a four-vector 



5 = ffoCTQ + g • (5p) = (5o,g)- 



(16) 



Rather than using the standard (ca,, ctj^, CTz) basis, we 
choose to rotate to ((T||, c^, CTz), the subscripts || and _L 
indicating respectively the directions parallel and per- 
pendicular to the internal field b. Defining the rotation 
matrix R(</') by 
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, (17) 



one has 

3 

5^, 51 ^w' ("^^f A'' ' (5^J = (5o,5|h5-L,5^) (18) 

3 

fj.'=0 

Expanding in harmonics - we also drop the four-vector 
indices 

K{ip, if') = Kf''' + cos{ip - (p')K^'''> + sm{if - (f')K^'''> + ..(20) 
In the above we have defined 

/O 
I 



ifi 

Vo 



2i^i 
iCa 
^2 
2i^i 
(21) 



and 



K(-) = 



/ 0- 

~K2 

iv:2 

V 



(22) 



For the purpose of calculating polarizations and spin cur- 
rents the higher harmonics play no role and are thus ig- 
nored. 

By using that gj}^ = ^gtq ~ 1 ^ ^5^) • cr and performing 
a rotation to the new spin basis, one can write Eq. ([5]) as 

dtg' = ^ i-Mg' + (No + Ni)(.g') + (N2 + N3){Kg')]+Ss. 

(23) 

The matrices appearing in Eq. (|23p read 



M = 



1 



r VF 





1 
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' 3rsf 
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(26) 



N, = — 



/ 1 
10 
10 

V 1 

(27) 
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where r* is the elastic quasi-particle hfe time, defined as 
11 1 

T* ~ T Tsf 



(28) 



which we now use for convenience of notation, but will 
be later incorporated into the proper transport time. Fi- 
nally, Sg is the source term due to the electric field. We 
take this to be along the x-direction, so that 



/ cos ifi \ 



Se = |e|i;F^^5e(2tanh(e/2r)) 



V 



■ COS(p — 

- sinG9— 







(29) 



Solving for the Sz spin current flowing along y, we ob- 
tain 
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J devpiPygz 



= / d. 



4 



4 
3rsf 



lU! 



2ma 
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liU 



2ma 



2ma 



i{pyg±) - {Pxg\\)) 
{gy) 

(30) 



i.e. the continuity equation result, Eq. under homoge- 
neous conditions. In the third line we have used Eq. (jl7p 
to set {gy) — {Pyg±) ~ {Pxg\\)- Similarly, one obtains the 
complete expression for the frequency dependent Sy spin 
polarization 



1 



3t, 



- iuj \ + 



1 



3t, 



2iu} 



-\ -1 



(31) 



Besides ^/Tsf, there appear in the above two other dif- 
ferent time scales 



Ttr T Tsf 



TE T Tsf 



The first, Ttr, is the total transport time. The second, 
te, is the generalization of the characteristic time related 
to the Sy spin polarization introduced in ([7]). By using 
Eq. ([3T|) in Eq. (PO]) . one obtains the expression for the 
frequency dependent spin Hall conductivity 



asH(uJ) = — 



Air V 3r, 



iuj ) 2{apF)'^ 
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FIG. 2: (Color online) Real part of the frequency depen- 
dent spin Hall conductivity in units of the universal value 
|e|/87r for appT = 1 (a) and appT = 5 (b). The different 
curves correspond to different values of the ratio t/ts/ = 
0,0.1,0.2,0.3,0.4 (from top to bottom at the maximum of 
Re asn)- 



Its real part is displayed in Fig. [5] for different values of 
the disorder parameter apFT. In the limit lu 0, the 
magnitude of the spin Hall conductivity depends on the 
value of apFT as well as on the ratio T/Tsf. In the ab- 
sence of magnetic impurities one has the known result 
o'sH = 0. As spin flip scattering grows, the conductiv- 
ity reaches values of the order of the "universal" |e|/87r. 
This was noted already in [l^ . where however, as pointed 
out in the beginning, angle dependent scattering was not 
considered. Large values of apFT can be achieved both 
in III-V and II- VI semiconducting materials. Doping the 
latter with Mn allows to control the spin-flip time Tsf 
while only weakly affecting the electrons mobilit y^^i'^^'^'^ , 
even though it is not perfectly clear whether these can 
appropriately be described in terms of the linear Rashba 
modelSl. Additionally, for certain frequencies one can 
see crossing points [cot « 0.5 and wr w 2 in Fig. [2] 
(a)] at which magnetic disorder has no effect on the spin 
Hall response. Such points are well defined only when 
apFT « 1. For clean {apFT ^ 1) or dirty {apFT 1) 
samples the different curves cross each other over a pro- 
gressively wider range of frequencies. 

Finally, in the diffusive regime, ujTtr ^ 1, apFTtr ^ 1, 
and assuming Ttr /Tsf ^ '^^TE/Tsf <ti 1, one obtains the 
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following spin-diffusion equations 



dtSx 



dtSy = 



dtSz 
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Sy — aNQ\e\S — 



(33) 
(34) 
(35) 



where {2apFTtr)'^ /2Ttr = I/t, is the Dyakonov-Perel spin 
relaxation rate, tied to Rashba spin-orbit coupling. From 
Eq. (|34p the sensitivity of the in-plane spin polariza- 
tion on spin-flip scattering is apparent: in the stationary 
limit the source (proportional to £) is balanced by the 
spin relaxation. Spin-flip scattering leaves the source un- 



changed, whereas it enhances the relaxation rate so that 
in the end Sy is reduced. 

In conclusion, we studied the combined effect of long- 
range and magnetic disorder on voltage induced spin 
polarizations and the related spin Hall currents in a 
Rashba 2DEG. We investigated homogeneous but non- 
static conditions, from the dirty {appT 1) to the clean 
{apFT » 1) regime. Care is required when treating long- 
range disorder because of the two-band structure of the 
problem, while magnetic impurities, even in low concen- 
trations, play a non-trivial role beyond that of a simple 
redefinition of the time scales. 
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